Classical search algorithm with resonances in y/N cycles 
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In this work we use the wave equation to obtain a classical analog of the quantum search 
algorithm and we verify that the essence of search algorithms resides in the establishment of 
resonances between the initial and the serched states. In particular we show that, within a set of 
N vibration modes, it is possible to excite the searched mode in a number of steps proportional to 

Vn. 
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I. INTRODUCTION 



It has been shown that a quantum search algorithm 
is able to locate a marked item from an unsorted list of 
N elements in a number of steps proportional to y/N, 
instead of proportional to N as is the case for the usual 
algorithms employed in classical computation. The most 
well studied quantum search algorithm is the one due to 
Grover Q , where the search is performed by alternatively 
shifting the phase of the searched for state, and amplify- 
ing its modulus. A continuous time version of the original 
Grover algorithm has been described by Farhi and Gut- 
mann 2]. In a recent work Q, we have presented an 
alternative continuous time quantum search algorithm, 
in which the search set is taken from the set of eigenvec- 
tors of a particular hamiltonian. The search is performed 
through the application of a perturbation which leads to 
the appearance of a resonance between the initial and 
the searched state. That study also provided a new in- 
sight on the connections between discrete and continuous 
time search algorithms. Our search algorithm can be im- 
plemented using any Hamiltonian with a discrete energy 
spectrum. However, we would like to emphasize that the 
possibility of establishing a resonance between two states 
is an intrinsic property of oscillatory motions in general 
and not an exclusive property of quantum mechanics as 
described through the Schrodinger equation, as in our 
algorithm. 

A link between quantum computation and classi- 
cal optical waves has been well established by several 
authors^ IE IE IE IE El- It is possible to simulate the 
behavior of some simple quantum computers using clas- 
sical optical waves. Although the necessary hardware 
grows exponentially with the number of qubits that are 
simulated, and thus these simulations are not efficient, 
optical simulations could still be very useful. In fact, 



as some simulations of quantum algorithms employ op- 
tical simulations, a classical analogue of such a quantum 
search algorithm might be a valuable tool to test the 
functioning of the optical system as a computer. 

In this work we present a continuous time search algo- 
rithm, which is controlled by a classical wave equation, 
showing explicitly, once again, that the search algorithm 
is essentially a resonancephenomenon between the initial 
and the searched states [E E3 • The paper is organized as 
follows: in the next section we briefly describe our quan- 
tum search algorithm. Then, in section 3, we develop 
the search model using the ordinary wave equation. In 
section 4 we present results obtained with this model. 
We end by discussing these results and extracting con- 
clusions, in the last section of this work. 



II. QUANTUM SEARCH ALGORITHM 

We consider a continuous time quantum search algo- 
rithm which is controlled by a time dependent Hamil- 
tonian H. The wavefunction \^{t)) satisfies the 
Schrodinger equation 



in 
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where H = H + V{t) and we have taken the Planck 
constant % = 1. Here Ho is a known nondegenerate 
time-independent Hamiltonian with a discrete energy 
spectrum {e n } and eigenstates {|n)}. V(t) is a time- 
dependent potential that shall be defined below. We 
then consider a subset N of {|n)} , formed by N states, 
which will be our search set. Let us call \s) the unknown 
searched state in N whose energy e s is given. We assume 
that it is the only state in N with that value of the energy. 
Therefore, knowing e s is, in our algorithm, equivalent to 
the action of "marking" the searched state, in Grover's 
algorithm. The potential V that produces the coupling 
between the initial and the searched states, is defined as 



3] 



V{t) = \p) (j\ exp (iuj sj t) + \j) (p\ exp (-iw aj t) (2) 



2 



where the eigenstate \j), with eigenvalue Sj, is the initial 
state of the system and it is chosen so that it does not 
belong to the subset N. Above, \p) = -7= \n) is an 

unitary vector which can be interpreted as the average of 
the set of vectors in N, and uj S j = ej —£ s - This definition 
ensures that the interaction potential is hermitian, that 
the transition probabilities W n j = | V^(t) |j) | = jj, 
from state \ j) to any state of the set N are all equal, and 
finally, that the sum of the transition probabilities ver- 
ifies ^ W n j = 1. The objective of the algorithm is to 



force with frequency Uj 



which may be derived from 



a potential V(t), the sought after pendulum s receives, 
after some time, most of the oscillation energy of the 
pendulum j, so it may be readily identified. This sys- 
tem, closely related to the one proposed by Grover [lOj, 
also has the property that the identification is performed 
through the gradual exchange of energy between the os- 
cillators, without a net external energy input. 

Returning to the quantum search problem, the wave- 
function appearing in Schrodinger's eq. may 
be expressed as an expansion on the eigenstates {\n}} of 
Ho, 



!*(*)) = ^ Om(t) exp (-ie m t) |m). 



(3) 
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FIG. 1: (a) Initial configuration of the system: all N pen- 
dulums at rest, except pendulum j. (b) Situation after y/~N 
cycles: most of the energy was tranferred to the searched 
pendulum s. 

find the eigenvector \s) whose transition energy from the 
initial state \j) is the Bohr frequency ou S j. At this point 
we consider important to make a digression to visualize 
the way the perturbation potential V(t) operates on the 
system. To better understand the process, let us consider 
a classical system composed of N pendulums of different 
lengths, and, consequently, different oscillating frequen- 
cies. The search problem, illustrated in fig. ^ consists 
of identifying the pendulum s with a given frequency of 
oscillation lu s . To perform the search, one suspends all 
the pendulums from the same beam, and puts to oscil- 
late one of the pendulums, j. By applying a periodical 



The time dependent coefficients a m (t) have initial condi- 
tions Oj(0) = 1, a m (0) — for all m ^ j. After solving 
the Schrodinger equation, the probability distribution re- 
sults in, 



Pj ~ cos 2 (ft t) , 
P s ~ sin 2 (ft t) , 
P n ~ 0, n ^ j and n =/= s, 



(4) 



where ft = 1/yN. From these equations it is clear that 
for t = 7r/(2ft) a measurement yields the searched state 
with a probability very close to one. This approach is 
valid as long as all the Bohr frequencies satisfy u> nm ^> ft 
and, in this case, our method behaves qualitatively like 
Grover 's. 



III. SEARCH ALGORITHM USING CLASSICAL 
WAVES 

In this section we build a classical search algorithm 
which is controlled by a classical wave equation with an 
added perturbation. The task of this algorithm is, start- 
ing from an initial oscillation mode, to excite the searched 
mode, whose frequency is known, in 0(y/~N) steps. The 
unperturbed wavefunction ty(x,t) satisfies the ordinary 
wave equation in one spacial dimension, 



d 2 ^{x,t) _ 2 9 2 *(a;,t) 



dt 2 



dx 2 



(5) 



where c is the constant wave velocity. This equation plays 
the same role as the Schrodinger equation in the algo- 
rithm shown in the previous section. From eq. J^J) we 
define the basis on which the algorithm is built. Let us 
call L the length of the region where the wave equation 
is solved. Then the canonical basis is given by the modes 
exp(— iujt) sm(Kx) where K is the wave number and u) is 
the circular frequency, connected through the dispersion 
relation c = uj/K. In this case, the wave number and 
the frequency are of the form K n — Kqh and u> n — ujQn 
respectively, where n is an integer, Kq = 2tt/L and u>q 
= cK . 
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As the analogue of the Schrodinger equation eq.Q) 
with the complete Hamiltonian H = H + V(t) we con- 
sider the wave equation with a perturbative force term, 



d 2 *(x,t) 2 d 2 ^(x,t) , d(V(t)T*(x,t)) 

~c ^-s h cluq (b) 



dt 2 



dx 2 



Ox 



where V(t) is a external time depend perturbative poten- 
tial. The solution of eq. © can be expressed as a linear 
combination of the oscillation modes, of eq. (JSJ 

*(x, t) = a n (t) exp(-iu> n t) sin(K n x), (7) 

n 

where the coefficients are 



0"n(t) = exp(io; n t)— / i) sin(_K„2;) dir. (8) 



In order for the algorithm to work, the operator T in 
eq.© is defined as 



T^!(x, t) = a„(i) exp(-iw„) cos(JsT ri 



(9) 



In order to obtain a greater connection with the pre- 
vious section, we shall now adopt the Dirac notation for 
the coordinate space |a;), and for momentum space \n). 
In this framework, the solution of the wave function can 
be expressed as 



I*) = ^ «n( f ) ex P j-iVnt) \n), 
n 

where we can write 



(10) 



a n (t) = 2ex P^"*) / ( /,./, ( |,.}/,.|*}. 



with (n\x) = sm(K n x). 

Once again, the task of the algorithm is to excite the 
oscillation mode |s) starting from the initial mode 
with the transition frequency u) s j = luj — uj s . To perform 
this task, we choose as perturbation potential to excite 
the resonance between the initial and searched states the 
classical analogue of the one given by eq. J2J ■ In our mixed 
notation, this time depended perturbation is 

V{t)#{x, t) = (x\V(t)\V) = £ {n\V(jt)\9) sin (K n x) , 



where 

(n\V(tm 



(12) 



{(1 - 8 n j) a,j(t) exp (iuj S jt) (13) 



S 7l j exp (—iwsjt) ^ ai(t) exp (— iwjt)}. 



The last equation has been obtained from eas. (|2l3() . Us- 
ing eas. 17112113(1 in eq.® and projecting the result on the 
mode n the following set of equations for the amplitudes 
a n (t) are obtained 



d 2 a n (t) 
dt 2 



da n (t) u> 
2iu n — — — = 7=Wn (1 - o n j) aj(t) exp (iuj sn t) 



N 



(14) 



+ WjS n j ai(t) exp (iuji s t)}. 



These equations will be solved in the next section but 
meanwhile we address the question of their qualitative 
behavior. The eas. (|14l> involve two time scales, a fast 
scale associated with the frequencies u) n , and a slow scale 
associated with the amplitudes a n (t). In a time interval 
larger than a characteristic time of the fast scale, the 
most important terms are those that have a very small 
phase, as the others average to zero. In this approxima- 
tion the previous set of equations becomes 



daj (t) 

dt 
da s (t) 

dt 
da n (t) 
dt 



—ifla s (t), 
—iClaj(t), 
for all n ^ s, j; 



(15) 



where Q = cjo/(2viV). These equations represent two 
coupled oscillators where the coupling is established be- 
tween the initial and the searched states. Solving them 
with initial conditions dj(0) — 1, a s (0) = we obtain 



aj (t) 
a s (t) 



cos(fit), 
i sin(fit). 



(16) 



It is important to note that the previous approach is 
valid only if all the frequencies verify u> n 3> The 
amplitude of the searched state is maximum at the time 
r = tt^/N /loq. This time is equivalent, in the Grover algo- 
rithm, to the optimal time to perform the measurement. 
Then, when the previous approximations are valid, our 
method behaves qualitatively like Grover's. The time r 
has a different prefactor than the searching time of the 
Grover algorithm, [l| but this difference is not important 
because the resonance potential eq.(J2Jl can be adjusted 
through multiplication by a constant factor. 



IV. NUMERICAL RESULTS 

We have integrated numerically eas. i|14|) with initial 
conditions aj(0) — 1, <x n (0) = for all n =/= j, and 
have verified that the solutions so obtained are indepen- 
dent of the remaining initial conditions for the derivatives 
daj ^ and da ;(°) , The calculations were performed using 



dt 



dt 
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N = 1 000 
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FIG. 2: Time evolution of the squared amplitude of the ini- 
tial state (dashed line) and the searched state (full line) for 
different number of modes N. The frequencies co n are taken 
with n = 1,2,3, ...,N. 



a standard fourth order Runge-Kutta algorithm. Choos- 
ing an arbitrary frequency for the searched mode, we fol- 
low the dynamics of the set N. We have verified for sev- 
eral values of N that the most important coupling is be- 
tween the initial and the searched mode; other couplings 
being negligible as discussed in the previous section. We 
worked with two types of sets N, in the first set the mode 
frequencies u> n = luqh are taken with n — 1, 2, 3, TV and 
in the second set with n = 1,4, 9, N 2 . The square of 
the amplitude of the searched and the initially loaded 
modes are shown for the first set in Fig|3 The tempo- 
ral evolution was normalized to the characteristic time 
r. Each panel shows the square of the amplitude for 
N = 50, 100, 1000 modes respectively. From this figure, 
it is clear that the amplitudes of the initial mode and 
the searched mode alternate in time, as in the Grover 
algorithm. 

The time at which the squared amplitude of the 
searched state is maximum and its value is very near 
one agrees with our theoretical prediction r. The agree- 
ment is much improved when N increases because the 
condition f2 << u) nm is better satisfied as f2 ~ 1/ \/~N. 
When t ~ 2r the amplitude of the initial mode shows 
strong oscillations which decrease as the inverse of the 
search set size. We also observe that the frequency of 



FIG. 3: Squared amplitude of the initial state (dashed 
line) and the searched state (full line) as a function of time 
for TV = 100 modes. The frequencies uj n are taken with 
n = 1,4, 9,. ..,7V 2 . 



these oscillations increases with TV. 

FigEl shows the squared amplitude of the initial and 
the searched modes as a function of time for the second 
set. In this case the calculation agrees with the theoreti- 
cal prediction for all N > 10, and thus we have displayed 
only the calculation for N = 100. 

The search algorithm dynamics for the two sets, are 
rather similar. The differences observed are due to the 
properties of their frequency spectra. While the fre- 
quency spectrum of the first set increases linearly with 
the mode number, the frequency of the second set has 
a quadratic growth. Then in the first set, an additional 
interference effect must be expected because many more 
frequencies u sn have the same value for different values 
of s and n in the case of the first set than in the case 
of the second set. In this last case the frequencies u> sn 
behave pseudorandomly and the additional interferences 
do not appear. 

V. CONCLUSIONS 

We have developed a search algorithm using a clas- 
sical wave equation with perturbation. This algorithm 
behaves like the Grover algorithm, in particular the op- 
timal search time is proportional to the square root of 
the size of the search set and the probability to find the 
searched state oscillates periodically in time. The effi- 
ciency of this algorithm depends on the spectral density 
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of frequencies of the set where the search is made. A 
larger separation between frequency levels maximizes the 
probability of the searched state and allows for a better 
precision. Nevertheless the classical implementation of 
this algorithm requires N oscillation modes, whereas in 
a quantu m sy stem the implementation will only need log 2 
N qubitsdfl. 

Presently most of the attention is devoted to the de- 
velopment of quantum computational devices. In partic- 
ular, there have been several recent attempts to employ 
optical simulations to run quantum algorithms ( 0,0,13 )• 
These studies suggest that such optical implementations 
are possible for small systems with the presently available 
technology. Here, we have developed a general classical 
analogue of a quantum search algorithm which operates 
in the same way as the one implemented in the quantum 



case through optical devices. Thus, our work not only 
helps to understand the essence of the quantum algo- 
rithm, but also may serve as a testing ground for practical 
implementation of optical computational systems. Cer- 
tainly much further experimental and theoretical work 
will be needed to develop the practical applications of 
this algorithm. 
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